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Abstract 

We describe electron transport through small metallic grains with Coulomb 
blockade effects beyond the perturbative regime. For this purpose we study 
the real-time evolution of the reduced density matrix of the system. In the 
first part of the paper we present a diagrammatic expansion for not too high 
junction conductance, h/Air'^e'^Rt ^ 1, in a basis of charge states. Quan- 
tum fluctuations renormalize system parameters and lead to finite lifetime 
broadening in the gate-voltage dependent differential conductance. We de- 
rive analytic results for the spectral density and the conductance in the limit 
where only two charge states play a role. In the second part of the paper we 
consider junctions with large conductance, h/Ae^Rt ^ 1. In this case contri- 
butions from all charge states, which broaden and overlap, become important. 
We analyze the problem in a quasiclassical approximation. The two comple- 
mentary approaches cover the essential features of electron tunneling for all 
parameters. 

I. INTRODUCTION 

Electron transport through mesoscopic grains is strongly influenced by the large charging 
energy, Eq = e^/2C, associated with the low capacitance C of the system ^-^. An inter- 
esting example is the "single-electron transistor" where a small metallic island is coupled 
via tunnel junctions to leads and via a capacitor to a gate voltage. At low temperatures, 
T -C Eq, a variety of single-electron phenomena have been observed in this system, including 
the Coulomb blockade and oscillations of the conductance as a function of a gate voltage. 
If the dimensionless conductance of the tunnel junctions between the island and the lead 
electrodes, 

— ^^ _ ^ (-\\ 

is small, on a scale defined by the quantum resistance Rk — 25.8 kQ, the charge of each 
island is a well-defined variable. In the limit at ^ 1, the sequential single-electron tunneling 
can be studied in perturbation theory [|T],^ ; and descriptions based on a master equation or 



equivalent simulations of the stochastic dynamics are sufficient to account for the dominant 
features observed in single-electron devices. 

Recent experiments beyond the perturbative regime show deviations from the classical 
description, e.g. a broadening of the conductance peaks much larger than temperature P,|l. 
This indicates that, in general, quantum fluctuations and higher-order coherent processes 
should be considered. Even in the limit of weak tunneling, at < 1, nontrivial features 
appear in the vicinity of the Coulomb blockade threshold, when two charge states become 
nearly degenerate and perturbation theory fails. Several theoretical papers IQ-Q dealt with 



the problem of higher-order processes, exploiting the physical picture of electron tunneling 
via discrete charge states. This includes "inelastic cotunneling" [0,|l3) where in a second- 



order process in the parameter at electrons tunnel via a virtual state of the island. An 
extension of this process, which gains importance near resonances, is "inelastic resonant 



tunneling" |]10|,|T3[, a process where electrons tunnel an arbitrary number of times between the 
reservoirs and the islands. The term "inelastic" indicates that with overwhelming probability 
different electron states are involved in the different steps of the higher order processes. The 
description can been extended to describe strong tunneling through single level quantum 
dots ^. 

If the conductance of tunnel junctions is not small, at ^ 1, the physical picture changes. 
In this case the inverse lifetime F = l/R^C and, hence, the broadening of the excited charge 
states due to quantum fluctuations exceed the typical level spacing of excited island states, 
hV ^ Eq. Thus charge levels overlap and the concept of tunneling via discrete charge 
states becomes ill-deflned, raising the question whether charging effects survive under such 
conditions or whether they are washed out completely by strong quantum fluctuations. In 



Refs. [|16| , p!7| ,p|, p!8| , p!9| , p!2| it was demonstrated that at sufficiently low temperatures even for 
large values of at quantum fluctuations of the charge do not destroy Coulomb blockade of 
tunneling, but they lead to a strong renormalization of the effective junction capacitance, 
Coff oc Cexp(27r^at). The exponential dependence on at had been derived independently 
by renormalization group arguments |]16|,|12|, instanton techniques [0, and Monte Carlo 



studies 1 12; 2^. One important consequence of the strong capacitance renormalization with 



increasing at ^ 1 is the exponential reduction of the temperature limit below which charging 
effects can be observed. 

This article is devoted to the calculation of the conductance of a SET transistor beyond 
perturbation theory in at, in a range of parameters which is accessible to experiments. 
The island contains a large number of electrons which are coupled strongly by Coulomb 
interactions. We, therefore, cannot proceed with ordinary perturbation theory. Rather, 
we reformulate the quantum mechanical many-body problem of these electrons in a real- 
time path-integral representation. In order to handle the Coulomb interaction we perform 
a Hubbard-Stratonovich transformation which introduces a phase as a collective variable. 
We trace out all microscopic degrees of freedom and arrive at an effective action of the 
system [0,^, similar in structure to that known from the studies of Ohmic dissipation in 
quantum mechanics |2^. This procedure is addressed in Section 2. 

After a change from the phase to a charge representation we are able to perform for 
at ^ 1 a diagrammatic expansion of the time evolution of the reduced density matrix. In 
a charge representation we can identify sequential, co- and resonant tunneling processes 
with certain classes of diagrams. A restriction to two charge states allows us to evaluate 



the spectral function and the conductance of the system analytically. The results will be 
presented in Section 3. At higher temperatures more charge states play a role, which in 
general requires a numerical study of the diagrammatic expansion. 

In the opposite limit of strong tunneling, at ^ 1, many charge states play a role, and a 
formulation in terms of the phase, which is canonically conjugated to the charge, is more 
convenient. This limit is discussed in Section 4. We analyze quantum dynamics of the 
phase variable in a semiclassical (saddle-point) approximation and obtain an expression for 
the system conductance valid at not too low temperatures T ^ e^/2Ceff. The exponential 
renormalization of the effective capacitance for strong tunneling widens this temperature 
range substantially. The two approaches cover the essential features of electron tunneling 
for all parameters. 

In Section 5 we review briefly some results obtained earlier within different imaginary 
time techniques, e.g. renormalization group and instanton methods, and compare these 
results with those of our real time analysis. 

II. FORMULATION OF THE PROBLEM 

We consider a metallic island coupled by two tunnel junctions (L,R) to two leads and 
capacitively to an external gate voltage Vg. An applied transport voltage V = Vi^ — Vr 
drives a current. A microscopic description of this single-electron transistor is based on the 
Hamiltonian, H = Hi^ + H^ + Hi + H^h + -f^t.L + -f^t,R- Here H^ = J2 ^karO-kar'^kar describes 

kcr 

noninteracting electrons in the left and right lead, r= L,R, and Hi = Yl ^qacl^Cqa^ models the 

go- 
island states. The Coulomb interaction is accounted for in a capacitance model 



Hch = Ec Ij^cl^Cq^-nA 



(2) 



The energy scale Eq = e^/(2C) of the transistor depends on the total island capacitance, 
C = Ci^ + Cn+Cg, determined by the left and right tunnel junction and the gate capacitance. 
The charging energy can be tuned continuously by the "gate charge" 

Qg = -erig = -(ClV^ + CrVr + C^Vg) . (3) 

The tunneling Hamiltonian Htr = J2 [T'^^cLk^^Cgfj + h.c.) describes tunneling between the 

kqa ^ ^ 

island and the left and right leads. The matrix elements are related to the tunnel conduc- 
tances by R;^ = {e'^/h)J2N^iO)N[{0)\T'''\'^, where A^(0) denotes the densities of states of 

the island and the leads, respectively. In the following we will consider "wide" metallic 
junctions with iV 3> 1 transverse channels. Extending the spin summation they can be 
labeled by the index a = 1, ...N. In the following we will put % = 1 (except when it enters 
the quantum of resistance). 

Our aim is to study the time-evolution of the density matrix. We shortly sketch the 
main steps of the derivation of this description: 

- The time evolution of the density matrix introduces two propagators, a forward and back- 
ward propagator, which get coupled when we trace out electron degrees of freedom of the 



reservoirs. The procedure is known from the work of Caldeira and Leggett ||22| who, gener 



ahzing earher work of Feynman and Vernon, studied the influence of Ohmic dissipation on a 
quantum system. Similarly the influence on electron tunneling was described in Refs. pT|,3. 
Here, we generalize the later work from a single tunnel junction to the transistor. 

- In order to describe the Coulomb interaction between electrons we introduce via a 
Hubbard-Stratonovich transformation the electric potential of the island V{t) as a macro- 
scopic field. The interaction between electrons is replaced in this way by an interaction with 
the collective variable. 

- We treat the leads as well as the electrons in the island as large equilibrium reservoirs. 
The electrochemical potentials of the reservoirs are fixed, fi,- = —cVt- for r = L,R. The only 
fluctuating field is voltage of the island V{t). The definition eV{t) = —ifit) relates V{t) to 
a phase ip{t). Its quantum mechanical conjugate is the number of excess electrons n{t) on 
the island. As a consequence of the procedure outlined so far, the macroscopic field n(t) is 
independent of the microscopic degrees of freedom described by Cg^j and cj^. At this stage, 
the electronic degrees of freedom can be traced out. 

- The time evolution of the reduced density matrix p(t; ipi, ip2), which depends only on the 
phase variable cp, can thus be expressed by a double path integral over the phases corre- 
sponding to the forward and backward propagators ipj (j = 1, 2) 

CO CO Ifili ip2f 

Pc{tf,(pii,ip2i) = / difii / d(p2i / V[ipi{t)] / V[ip2{t)]exp{iS[(pi{t),(p2{t)])pc{tu(piuip2i) ■ 

— oo — oo ipu ip2i 

(4) 

- The form (^ describes the situation where charges can take any continuous value and the 
phase is an extended variable. However, in our physical system the charge on the island is 
quantized in units of the electron charge e. In this case the phase variable is compact (i.e., 
the states ip and ip + 27i are equivalent), and we rewrite @, introducing integer winding 
numbers mi,m2 = 0, ±1, ±2, . . ., 

oo oo ipi{+2nmi ip2[+2nm2 

Pd{tf,'Pli,f2f) = Y. / '^V^li / '^'^2i / ^blW] / 1^[f2it)] 

exp {iS[(p>i{t), ip2it)]) pd{tu v^ii, V22i) • (5) 

The two integrations can be combined to a single integral along the Keldysh contour, 
which runs forward and backward between ti and tf along the real-time axis. As a result the 
reduced propagator H is written as a single path integral along this contour 



H = tr 



Po^KCxp i-i dtH{t) 



I?[^(t)]exp(z5[^(t)]). (6) 



Here the collective variable (f{t) and the time integral are defined on the Keldysh contour 
K, and the time-ordering operator Tk orders the following operators accordingly. 

The effective action entering the propagator is ^[(^^(t)] = Sch['p(t)] + Stlipit)]. The first 
term represents the charging energy 



SchHt)] = / dt 



c (m 



+ %0(^) 



(7) 



Electron tunneling is described by St[(p{t)\, which, in the case of wide metallic junctions, is 
expressed by the simplest electron loop connecting two times. 



StW)] 



2vriE 

r=L,R 



dt / dt'a^{t,t')e 

K JK 



l\fJ-f{t)Q-iv{t') 



(8) 



The kernels a^(t,t') = (y.f{t — t') for t < t' (t > t') depend on the order of the times along 
the Keldysh contour. Their Fourier transforms are |^,10,13 



af{u) = ±at,r 



LU 



/ij. 



exp[±(cj - iJ,j.)/T] - 1 



(9) 



They are proportional to the dimensionless tunneling conductance at,r = h / {At^"^ e^ R^) be- 
tween the island and the leads r = L,R. 

For large systems, the phase behaves almost like a classical variable while its conjugate 
variable, the charge, fluctuates strongly. A natural basis is then the phase representation. 
In the presence of strong Coulomb interaction, however, the situation is different: the phase 
underlies strong fluctuations while the time evolution of the charge is almost governed by 
classical rates. For this reason, it may be useful to change from the phase to the charge 
representation. The time evolution of the density matrix in a charge representation depends 
on the propagator from ni forward to n'^ and on the backward branch from n'g backward to 



n2- It is given by the matrix element of the reduced propagator |13 



n2,n2 



d^i / d^\ / d^'^ / ci<^2e^"^^^e-™i^ie^'^2^2e-*"2^2 
1 



(10) 



1/32,1/32 



In the charge representation the charging energy is simply described by Sch[n{t)] 



III. EXPANSION IN THE TUNNELING CONDUCTANCE 



A diagrammatic description is obtained by expanding the tunneling term exp (i5't[<^(t)]) 
in the reduced propagator and integrating over (p. Each of the exponentials exp[±z(/9(t)] 
describes tunneling of an electron at time t. These changes occur in pairs in each junction, 
r=L,R, and are connected by tunneling lines Q^{t,t'). Each term of the expansion can be 
visualized by a diagram. Several examples are displayed in Fig. |I]. The value of a diagram is 
calculated according the rules which follow from the expansion of Eq. (|T^) and are presented 
in detail in Ref. [|l^ . 



n! 



The propagator from a diagonal state n to another diagonal state n' is denoted by 

n' 

'„/ = Unn'- It is the sum of all diagrams with the given states at the ends and can 



be expressed by an irreducible self-energy part E„„', defined as the sum of all diagrams 
in which any vertical line cutting through them crosses at least one tunneling line. The 
propagator can be expressed as an iteration in the style of a Dyson equation, n„„/ = 

^n^Snn' + Y.'^n n" ^n" n' '^n' ■ The term H^^^ describcs a propagation in a diagonal state 

n" 

which does not contain a tunneling line. The stationary probability for state n follows from 
Pn = J2 Pfi' n„/_„ (in which P^'> is the initial distribution) and is not the equilibrium one if 

n' 

a bias voltage is applied. Our diagram rules then yield 

= 2_,[~-PnSn,n' + P-n'^n' ,n] ■ (H) 

n' 

We recover the structure of a stationary master equation with transition rates given by S„/^„. 
In general, the irreducible self-energy E yields the rate of all possible correlated tunneling 
processes. We reproduce the well-known single-electron tunneling rates by evaluating all 
diagrams which contain no overlapping tunneling lines. Similarly cotunneling is described 
by the diagrams where two tunneling lines overlapping in time, as shown in Fig. |T]. 

We calculate the current I^ fiowing into reservoir r = L, R by adding a source term to the 
Hamiltonian and then taking the functional derivative of the reduced propagator with respect 
to the source. The result Jr = —iejduj{a^{uj)C^{uj) + a^{ijj)C'^{ijj)} is expressed by the 
correlation functions C>(t,t') = -i(e-^^(*)e*^(*')) and C<(t,t') = i(e^^(*')e-^'^(*)) describing 
charge transfer at different times. These are related to the spectral density for charge 
excitations on the island by 2T{iA{ijj) = C^{u) — C^{uj). 

For sequential tunneling, the current reduces to 



with 



fi J ^, a[LU) 



A^'\^) = Il[Pn + Pn+i]S{u; - AO (13) 



and An = Eci,{n + 1) - E^hin) = Ec[l + 2{n - n^)]. 

At the minima of the Coulomb oscillations the system is in the Coulomb blockade regime, 
and cotunneling processes determine the conductance. The second order processes are de- 
scribed by diagrams as shown in Fig. |I[ A careful analysis of our diagrammatic expansion 
not only reproduces the known limits but also provides the needed regularization of di- 
vergences. We, furthermore, obtain new terms which are essential at the resonance. These 
results will be presented in a forthcoming publication ||14|| . 

At the resonance we have to include processes of arbitrary high order, since the process of 
resonant tunneling is essential. For definiteness, we concentrate on situations where only two 
charge states, n = 0,1, need to be considered. This is the case when the energy difference of 
the two states Aq = -Ech(l) — -Sch(O), the bias voltage eV = eVi^ — cVr, and the temperature 
T are low compared to Eq- If, furthermore, we restrict ourselves to matrix elements of the 
density matrix which are at most two- fold off-diagonal []13[, we can evaluate - in a conserving 



approximation - the irreducible self-energy analytically. The following results are derived in 
this limit. 

6 



Using the notations a^^uj) = a+(u;) + a^. (u) and a^u) = J2'^r{^), we find Pq = A_ and 

r 

Pi = A+ with A± = J duj a'^{uj)\n{uj)\'^ and 

7rH = [c.-Ao-aM]-i , aico) = f du;' ^^^yl— . (14) 

J UJ — UJ' + 20+ 

Again, the current is given by Eq. (^), but the spectral density becomes 

Af \ — (y{uj) . . 

^""^ ~ [cu-Ao-Rea(cu)]2+[Ima(cu)]2 " ^ ^ 

The following results depend on the parameter 

which also defines the parallel tunneling conductance 1/Rt = J2 V-Rr- In lowest order in 

r 

at we have A^^\ij) = 5{lu — Aq), and the classical result is recovered. In general, quantum 
fluctuations yield energy renormalization and broadening effects, which enter in the spectral 
density via the complex self-energy cr{u) given in Eq. (0). In order to evaluate (j{uj) we 
introduce a Lorentzian cut-off which we choose equal to Eq (since the energy difference to 
charge states which are not taken into account here is of the order of the charging energy). 
In this case we find 

Rea(^) = -5:at,.(^-/.r)[21n(^)-2RevI/(2^^)] (17) 

and lm.a{uj) = — 7ra(u;). The effect of the quantum fluctuations can be estimated from the 
spectral density in the limits T ^ eV, \uj\ or eV S> T, |c<j|. Then, the spectral density is 



^M= ,.. ...T:. ..... . (18) 



Z'^aico 
[io - ZAo]2 + [irZaiuj)] 

with 



Z-^ = 1 + 2atln{Ec/ msix{eV/2,2'KT}). (19) 

We observe a renormalization of Aq and at by Z and a broadening given by 7rZa{uj). 
From this result we conclude that lowest order perturbation theory is sufficient for 
atln{Ec/ Taa.x{eV/2,2nT}) <^ 1. At larger values, our results for resonant tunneling show 
clear deviations from sequential tunneling. 

A pronounced signature of quantum fluctuations is contained in the differential conduc- 
tance G = dl/dV. In Figs. | and | we present our results for the differential conductance 
in the linear response regime {V = 0). They clearly display the effect of resonant tunneling: 
- For comparison, we show on the left hand side of Fig. ^ plots which are obtained from the 
master equation description of sequential tunneling. 



G{T, n, 

Cfas 



E„exp 



T 



— ^[n — n„ 



oo 

E 



exp 



' T 



n 



n„ 



f^(l + 2(n-ng)) 



exp 



T 



(l + 2(n-ng)) 



(20) 

The asymptotic high-temperature conductance is Gas = ^/{Rl + Rr}- At low temperatures, 
when processes involving only two charge states dominate, the maximal classical conductance 
saturates at one half of the asymptotic conductances at high temperatures. The width of 
the peaks scale linearly with temperature. 

- The situation changes when resonant tunneling processes are taken into account (see the 
plots on the right hand side of Fig. ^. The maximal conductance and the peak width are 
renormalized by Z and Z^^ which depend logarithmically on temperature. For this reason, 
the conductance peak does not reach one half of the high temperature limit and decreases 
with lower temperatures, while the peak width is increased compared to the lowest order 
perturbation theory result. For an estimate of the maximal conductance, we use can the 
spectral density in the form of Eq. (|THD and perform the integral Eq. (|T^) analytically, 

-^-'y-^)] , (21) 

(The results shown in Fig. (^), however, were obtained by numerical analysis based on 
Eq. (pD.) 

Recent experiments [^,^ in systems with junctions with small barriers show, indeed, a 
broadening and decreasing height of the linear conductance peaks, which cannot be explained 
by thermal smearing and qualitatively agrees with our theory. 

The effects of quantum fluctuations are even more pronounced in the nonlinear differ- 
ential conductance when the transport voltage dominates over temperature. In Fig. ^ we 
compare the results of perturbation theory and resonant tunneling at T = assuming that 
for eV < 2Eq only two charge states n = 0, 1 are involved. 

- The sequential tunneling result for a symmetric transistor with at,L = c^t.R and Cl = Cr 




IS 



Gas 



E2(i_2n,)2 + 



[eVf/A Q f eV 



1 



n„ — 



[eVf ^ \AEc '^^ 2 J ' ^^^^ 

As a function of Ug it shows a series of structures of width CV/e with vertical steps at its 
edges. The width scales linearly with bias voltage. 

- Resonant tunneling leads to a renormalization of the height and width by Z and Z~^ 
respectively, which depends now logarithmically on the voltage (see Fig. ^. For this reason, 
the height of the structure is below the sequential tunneling result and further decreases at 
lower voltages, while the width is enhanced. Furthermore, the sharp edges are smeared out 
even in the absence of thermal fluctuations (since T = 0). 



IV. STRONG TUNNELING 

If the junction conductance is high and hence the fluctuations in the charge are strong 
the phase representation outlined above is a more suitable starting point for the analysis of 
the problem. It turns out that the dimensionless conductance appears in the form 



8 



which differs from the expansion parameter at of the weak tunneling expansion by a factor vr^. 

The real-time path-integral technique discussed above provides an expression for the reduced 

density matrix p{ipi,ip2)- If the island charge can vary continuously the density matrix 

+00 
Pc{<^i, ^2) is given by Eq. (^. It obeys the standard normalization condition / dippc{(p, (f) = 

—00 

1 and its time evolution is governed by the action S[(p{t)] = Jj^ dt^ i'^) + ^t[fit)]. In this 
case the charging energy does not depend on the gate charge. In the SET transistor, another 
physical situation is realized, where the island charge is discrete and quantized in units of e. 
This situation is described by the density matrix Pdi'^i, ^2)^ Eq- (ID, with a compact phase 

variable [0. The normalization of the density matrix (|^) is given by / d(fpd{f,f) = 1- It is 

— 7r 

sensitive to the gate charge. The comparison of Eq. (^) and (^) shows the following relation 

Pd(<^, ^) = A-' J2 e*'™^(™i-™^ Vclv' + 27rmi, cp + 27rm2) . (24) 

mi,r?i2 

+00 ^^ 

Here A is a normalization factor A = Y, I dipe^'^'^"'^"^ pd^p + 27rr7i, ip). The relation (^) can 

m -00 

also be used to establish the connection between the expectation values of physical quantities 
for systems with discrete and continuous charge distributions. The expectation value of an 
operator O{0), of the discrete-charge system, which is 27r-periodic in cp is 

n +00 

6)^ = /rf^O(v^)pd(^,^) = E /rf^O(^)e^2-"^™pc(^-27rm,^) 

E (0(^)e*2'^("s-")'") 
= ULl 11. (25) 

m 

Here we used an obvious identity J dip{...)pc{'p — 27im,(p) = ((...)e~^'^"^'^)c. 

Now we are ready to evaluate the tunneling current through a SET transistor. We first 
derive an expression for the expectation value of the current and then evaluate it with the 
aid of Eq. (]25|). The first part of this program will be carried out within the quasiclassical 
Langevin equation approach ||2l|j23H 26| derived under the assumption that fluctuations of 



the phase variable are weak. This assumption is justified if the fluctuations of the charges 
are strong. 

In a semiclassical approximation we find for the current through the the left and the 
right junctions r = L, R (see Refs. [21,^| for further details) 



^r = cA + ^^-U^). (26) 

e itr e 

It depends on the fluctuating voltage differences across the junctions, 0r/e = U — V^ for 
r = 1,2. Here U is the voltage of the island. Each current is the sum of displacement 

9 



current on the capacitor, a deterministic tunneling current and a shot noise contribution. 



The latter can be expressed as a state-dependent noise in the form ||21 



where ^r 



6 = 61 (t) cos v^r + 62 (t) sin Lf, , (27) 

1,2) are two independent Gaussian stochastic variables with correlators 



(^ri(0)6'j(t)) = 5r,r"^ij 



1 



duj 
— ( 

27r 



!^'^*CJ coth 



2T 



(28) 



We consider situations where the external impedance is negligible. In this case the 
phases ipr are linked to the transport voltage V hy ipi, — (f-n = eV. We can further assume a 
symmetric bias Vl = —Vr = V/2. In this case the voltage on the island can be expressed as 
U = (0L + 0R)/2e = Vg + ipg/e. Here 0g/e is the voltage across the gate capacitance. Finally, 
charge conservation implies ^l + ^r = —CgCpgje. Combining these relations with ( P^D we 
arrive, after averaging over the stochastic variables ^, at the expression for the expectation 
value of the current 



(^l) 



Rl + Rr 



(29) 



This expression will be evaluated further with the aid of relation (]25|) . If the fluctuations of 
the charge can be treated as Gaussian the contribution of the m-th winding number to the 
expectation value (^Sf ) can be estimated as 



i2nhm\ 
e ) ~ 

/ c 



exp (— 27r^ ((5r?,^ ) m^ 



(30) 



Thus provided that the charge fluctuations are not small {Sn"^)^ ^ 1 it is sufficient to retain 
in the expression (|25|) only terms with winding numbers m = 0, ±1. In this approximation 
we obtain IE5i 



^r ) +2 Ur cos(27r(n — rig)) 
l + 2(cos(27r(n-?2g)))^ 



(31) 



In the quasiclassical limit considered here the further analysis requires standard noise 



averaging of the solutions of Eqs. (^61) . As these equations are nonlinear in the phase, 
the exact solution cannot be constructed in general. In a wide parameter range, however, 
it is sufficient to proceed perturbatively in the noise terms. Making use of (p9D, (pl|) and 
assuming the phase fluctuations to be small \Sip\ ^ tt we arrive at the following expression 
for the current 



I{V) = G^V - loiV) 

-F{T,V) 



— e 



[h{V) - 2Io{V)\ cos + hiV) sm 



(32) 



Here 
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Qav = (e(n - rig)). 



Cl-Rl — C'r-Rr 



V + aK 



-Rl + -Rr 
is the average charge of the metalhc island. We further introduced the integrals 



loiV) 



2e 



7r(i?L + Rr) 



dt 



f ttT 



V sinh nTt 



e ^ ' 'K[t) cos I „,„ „ , 1 sm 



2(/?R + /?l) 



(33) 



(34) 



/i(\/) ^ 4e 

/2(V^) 7r(i?L + i?R) 



dt 



/ ttT 



VsinhTrTt 



-iy(i,y) 



^Wlt^-^^'^)) 



+A^Ck{tf^^^^{u{t,V)) 



cos fe{RR-RL)Vt\ . /e\/t 

sm 



sin \^ 2(i?R + i?L) ; V 2 
and we defined K{t) = i?t^(t)(l - exp(-t//?iC)), 1/Rt = 1/R^ + 1//?r and 



(35) 






vrT 



sinh7rT(ti — ^2) 



K{tut)K{t2,t) 



V — cos (^^^^^^^^^^-^ 
-Rr V -Rl + -Rr 



r=L,R 



(36) 
(|3|) should 



with K{t', t) = K {t') — K {t' — t) . A principal value of the time integrals in ([ 
be taken where needed. 

The function F{T, V) = 27i'^{6n'^) determines the temperature and voltage dependence of 
the charge fluctuations in the Gaussian approximation. It is given by an expression similar to 
W{t,V) (0) with the substitution K{ti,t)K{t2,t) -^ {ATT^C^/e^)k{ti)k{t2). Analogously 
u{t,V) is defined by (|3|) after the substitution K{ti,t) -> -{AnC^/e^)k{ti). 

We can simplify these expressions observing that in the limit of sufficiently high temper- 
atures and/or voltages 



maxleV, T} > wq = ^ — exp(— 2at + 7), 



TT^ 



(37) 



(here 7 = 0.5772... is Euler's constant) the results can be simplified further, since the time 
integration in ( p6D is effectively cut off at short times. Since W{t = min{l/T, 1/eV}, V) <^ 1 
we can set in the parameter range (|37|) in leading order approximation W{t,V) = 0. With 
this simplification the above integrals can be evaluated analytically. We obtain 



Io{V) 



e-Ro 



-Rl + -Rr 



E ■ 

r=L,R 



eRrV 



7r(i?L + i?R) 



Re^ 1 



1 



eRrV 



— ^ y Im^ I 1 



2ttTRoC 27rT(i?L + Rr] 

T. . ( eRrV \ 

-Re^ 1 - i 

1 



27tT{R^ + R^)J 

eR,V 

2-kTRqC " ^27rT(i?L + Rr) ^ 



(38) 
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and 

f(T.V) = F(0,0) + ?^ + 4^1„^ ' 



^ r=L,R 



e2 emo \27rTRoCj 

1 pRV M 

• (39) 



1 r^_,^M^y\^(,^ 1 . ei?.i/ 



i?L V Rl + Rr J V 2ttTRoC 2ttT{Rl + Rr) 



Here \E'(x) = r'(x)/r(x) is the digainina function. 

The last expression determines the temperature and vohage dependence of the charge 
fluctuations (5n^) at not too low T and/or V. At T = and V = the integral over time in 
F[T, V) diverges logarithmically at high frequencies. This divergence indicates a failure of 
the quasiclassical Langevin equation in this limit. The problem can be cured by observing 
that in thermodynamic equilibrium (zero voltage) in Gaussian approximation 

e-^(°'°) = (cos(27rn)) = ^ d^P;^(^n + ^^^) 



The expectation value involving the equilibrium density matrix can be evaluated |2S| with 
the result 

F(0,0)~2af (41) 

The functions Ii{V) and hiV) cannot be evaluated analytically even in the limit (p7D. 



Due to a fast decay of the exponential factor exp[— F(T, V)] in (|32D with increasing V and 



T it is sufficient to evaluate /i and I2 in the low voltage and temperature limit. In this limit 



the integral (|35|) reduces to 



1 22 
hiV)-2Io{V)^gG^V, g=^^ + 11.29. (42) 

at 

whereas the function l2{V) turns out to be small l2{V) ~ l^^ ~ and will be neglected 
below. 

We thus arrive at the following result for the I-V characteristics of a SET transistor 



I{V) = G^V - I„{V) - gG^l^e-"!'"-"! cos ( 



2^Q,An] . (43) 



The current is reduced below the classical result Gas^ by an amount loiV) and is modulated 
in a periodic way by the gate voltage. In the limit considered the modulation is a pure 
cos-modulation. The result (^) also describes the oscillatory behavior of the current as a 
function of the transport voltage, which is usually referred to as a "Coulomb staircase" . The 
amplitude of these oscillations decays exponentially with increasing voltage and temperature. 
We also recover the fact that the Coulomb staircase is pronounced only in asymmetric SET 
transistors. In a symmetric case the transport voltage drops out from the expression for the 
gate charge (^ . The I-V characteristics (|43|) is depicted in Fig. ^ for different temperatures 
and values of the gate charge. 

The linear conductance of a SET transistor can be easily derived from Eq. (^3]) in the 
hmit V ^0. We find 
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^21 = 1 _ /(T) - ^e-^(^'°) cos (27™g 



where n^ = CgVg/e and 
fiT) = ^ 



2atEc^'/, 2atEc\ , /, 2at^c 
7 + — ^^ 1 + — ^ +^ 1 ' 



TT^T V n^T J V TT^T 



(44) 



(45) 



These results are displayed in Figs. ^ and |^ in the temperature range T > IOwq, where we 
estimate the approximations used above to be justified. 

In the high-temperature limit the conductance becomes independent of the gate charge, 
but due to charging effects it is still reduced below the asymptotic value by 

For high temperatures this expression is valid for all (including small) values of at- The first 
nontrivial term in this expansion does not depend on at- The coefficient of the square term 
also contains an at -independent contribution k. Within the approximation W = used 
here we have k = 0. An improved approximation is obtained by expanding in W{t, V), or 
alternatively by treating the general expression for the system conductance perturbatively 
in at and then expanding in Eq/T. This procedure yields k, = 1/15, which for large at can 
be neglected against the first term 6({3)at/n'^. 

At lower temperatures the conductance is further suppressed by charging effects and it 
can be modulated by the gate charge Qg. In the figures the minimum and maximum conduc- 
tance values are presented corresponding to Qg = and Qg = e/2, as well as the Qg-averaged 
conductance. The modulation with Qg becomes more pronounced as the temperature is low- 
ered, however, it is exponentially suppressed with increasing dt (cf. Figs. ^ and |^. For 
at ~ 4 the modulation effect can hardly be resolved while the overall suppression of the 
system conductance G is very pronounced. 

Although the validity of the Langevin description description is restricted to high temper- 
atures and/or voltages, T, eV ^ wq, the validity range rapidly expands as dt increases. E.g. 
for the parameters Eq ~ 1 K and at ~ 2, we get wq in the range 15 niK. Further increase 
of at rapidly brings Wq below 1 mK. Therefore we can conclude that in the strong tunneling 
regime at > 2 -t- 3 our theory covers the experimentally accessible temperatures. Indeed a 
quantitative agreement without fitting parameters exists between our results ( ^^ - |46|) and 
those of the Saclay group [^] in the high temperature regime. For lower temperatures the 
quasiclassical Langevin equation approach can be applied only to sample 4 of Ref. f^ with 
at — 1.8. Other samples studied in Ref. have substantially lower conductance, and their 
low-temperature behavior should be described by the expansion in at presented in Section 
3. 

V. DISCUSSION 



In a number of earlier papers [|T^,0,|1^,|T^ the combination of charging and strong tun- 



neling effects in metallic junctions has been analyzed within imaginary time approaches. In 

13 



the limit of strong tunneling, at ^ 1, a. renormalization group equation for at can be derived 



dat/dlnuc = /5(at) , (47) 

where in the lowest order in at one has (3{at) = 1/2. Already this scaling approach captures 
the tendency of the effective junction conductance to decrease with decreasing T due to 
charging effects. In order to see that one should proceed with scaling from ujc ~ Eq to 
ujc ^ T and identify the (dimensionless) junction conductance with the renormalized value 
ati^c ~ T). This approach is sufficient for strong tunneling at high temperatures, namely 
if the final renormalized tunneling conductance still satisfies at{u}c ~ T) ^ 1. In general 
the strong tunneling approach may lead to a small renormalized conductance such that ( ^Tf ) 
ceases to be valid. For weak tunneling other scaling approaches, derived in an expansion 
in the tunneling conductance and equivalent to what we described in Section III, can be 
applied. In this situation, Falci et al. [jl2| suggested a 2-stage scaling procedure, where 
the renormalized conductance after the strong tunneling rescaling was used as an entry 
parameter for the weak tunneling scaling. 

Various theoretical approaches led to the conclusion the strong electron tunneling at ^ 1 
reduces the charging energy, i.e. the effective capacitance is renormalized. Panyukov and 
Zaikin [^ treated the problem by means of instanton techniques. They concluded that 



electron tunneling affects both the scale and the functional dependence of the ground state 
energy E{Qg). At not too low temperatures T ^ w they find 

E(Q,) = -^cos(^ (48) 



with 18 



w 



2>2atEc 



7]-2 



exp(-2at + 7)- (49) 



A similar result, differing only in the numerical coefficient, has been obtained in a semi clas- 



sical analysis of the effective action ||T2 



At lower temperatures the form of the lowest energy band E{Qg) turns out to be even 
more complicated 0,113] and the at- dependence of the prefactor of the expression for w 



changes from linear in at for T > w to quadratic in at for T = 0. Instanton techniques ||T8[ 
yield 

£^c,eflf oc a^Ecexp{-2at) . (50) 

The exponential dependence on at has been confirmed by renormalization group argu- 
ments |lT6|,p|,p!2| as well as Monte Carlo methods [0,0 . The prefactor remains a point 
of controversial discussions in the literature PDI. Irrespective of this detail an important 



consequence of the strong capacitance renormalization for at ^ 1 is the exponential reduc- 
tion of the temperature range where charging effects are observable. 

With the aid of relations (^), (^9|) we can derive the first order correction in 1/at in the 
renormalization group equation (|4^) |^ 
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/3(«t) = 1/2 + l/4«t. (51) 

This result has also been derived by direct RG methods [p8|] . 

A consequence of the renormalization group approach ( ^Tf ) has been pointed out in Ref. 
||18|| . It relies on the assumption that the system linear conductance is determined by the 
renormalized value Q;t(ci;c ~ T) as 

G'=^at(^c~T). (52) 



Combining the above scaling approach, the high temperature expansion ( ^Bj ) (with k = 0), 
and the expression ( ^1]) for /? to first order in I /at we get for the Qg-averaged conductance 

G If/ Aaf EA r 1 / Aaf ^cM 1 . ^ 

-— = 1 hnil + — '- — -— +nl + n 1 + ^ — r— }. 53 

Gas 2at\ \ 3il + 2&t)T) [ 2a, \ ?>{l + 2a,)T)\\ ^ ' 

Although the above scaling approach to the conductance calculation is intuitively attractive 
(and the result ( |5B| ) fits reasonably with the available experimental data |^,^) it has to be 
stressed that it depends on the unproven assumption (|52D . 

In contrast, the real-time path-integral techniques presented here are free from this ambi- 
guity and allow for a direct evaluation of the I-V characteristics and the system conductance. 
We note, furthermore, that the results obtained within the real and imaginary time methods 
are consistent with each other. E.g. the renormalization of the effective energy difference 



between the two lowest charge states, derived in Ref. |]12[ , is contained in the self-consistent 
solution presented in Section 3. Furthermore, comparing the expressions for wq ( P7| ) and the 
bandwidth w ( ^Qf ) we immediately see that these two parameters coincide up to a numerical 
coefficient: w = IGwq. This means the requirement for the validity of the quasiclassical 



Langevin equation (|37D roughly coincides with the requirement that the temperature (or 
voltage) is larger than the effective bandwidth w. 

Still no quantitative theory for the conductance at lower temperatures and not too low 
values ttt ~ 1 has been provided. Although the two limiting descriptions presented here do 
not allow for a quantitative description of this parameter range it satisfactory to notice that 
both show the same qualitative trend in this range. 

Another question of interest is the conductance at very large at S> 1 and very low 
T ^ Wq. In the limit at ^ 1 the conductance oscillations with Qg are exponentially small 
(cf. (H)). Then for all Qg from (||,|I|) we have 

G{T^wo)/Gas-b/at, b ^ 1. (54) 

Thus we can conjecture that the low temperature maximum conductance of a SET transis- 
tor is universal in the limit of large at being of the order of the inverse quantum resistance 
unit 2e^/7r/i. This conjecture is also consistent with the scaling analysis of Refs. |T6| , p!8|jr^ 
combined with the results of Section 3. Starting from large at ^ 1 we first use the renor- 
malization group procedure (|i7|J^ ) which should be cut at at(^c) ~ 1. In the second stage 



we expand in at ~ l/vr^ as described in Section 3 - starting with the renormalized value 
instead of the bare one. Apart from logarithmic corrections we thus arrive at the maximum 
conductance of order of the inverse quantum resistance, no matter how large the initial 
conductance is. 
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VI. CONCLUSIONS 

In this paper we have described single-electron tunneling in systems with strong charging 
effects beyond perturbation theory in the tunneling conductance. For this purpose we con- 
sidered the real-time evolution of the reduced density matrix of the system. We presented 
two approximation schemes: 

In the first part, valid for not too strong tunneling, at ^ 1, we presented a systematic 
diagrammatic expansion, which allowed us to identify the different contributions, sequen- 
tial tunneling, inelastic cotunneling and inelastic resonant tunneling. When we restricted 
ourselves to diagrams corresponding to maximally two-fold off-diagonal matrix elements of 
the density matrix we can formulate a self-consistent resummation of diagrams. At low 
temperatures we, furthermore, can restrict our attention to two consecutive charge states. 
In this limit, there exist no crossing diagrams, and we can evaluate the summation in closed 
form. The most important results are a renormalization of system parameters and a life-time 
broadening of the conductance peaks. These two approximations are justified for tunneling 
conductances satisfying at In [Eq/ max{e\//2, 27rT}) < 1 and allow for a qualitative analysis 
of the system conductance also for larger values of at- 

In the second part of the paper we developed an alternative approach based on quasi- 
classical Langevin equations for the junction phase ip. This approach assumes that fluctu- 
ations of the phase are small and that the noise can be treated perturbatively. This is a 
suitable approximation for large values at = T^'^ott or in the high temperature limit. For 
weak tunneling at ^ 1 this scheme turns out to be justified only for high temperatures 
and/or voltages max(T, eV^) ^ Eq,^ whereas for stronger tunneling, ctt ~ 1, phase fluctua- 
tions are substantially suppressed. The results derived in this approach are valid, provided 
max{T, eV} ^ at-Ec6xp(— 2Q;t)- This range expands rapidly with increasing otf 

In conclusion, we found an effective action description of a single-electron transistor. 
We analyzed it in two limits. The charge representation, which is valid as long as at ~ 1, 
provides the basis for a systematic diagrammatic description of coherent tunneling processes 
including resonant tunneling. The phase representation is suitable at large values of at ^ 1. 
In both cases we calculated the gate-voltage and temperature-dependent conductance of a 
single electron transistor. The dimensionless parameters in the two limits differ by a factor 
vr^at = at. As a result the range of validity of the two approaches overlaps and, at least 
qualitatively, the two approaches cover the whole range of parameters. 

The authors are grateful to D. Esteve, G. Falci and G.T. Zimanyi for useful discussions. 
We thank the members of the Saclay group for sending us their data prior to publication. 
The project was supported by the DFG within the research program of the Sonderforschung- 
bereich 195 and by INTAS-RFBR Grant No. 95-1305. 
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FIGURES 
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FIG. 1. A diagram showing from left to right: sequential tunneling in the left and right junction, 

a term preserving the norm, a cotunneling process, and resonant tunneling. 




FIG. 2. The linear differential conductance normalized to the high temperature limit. 
Left hand side: result from a master equation in lowest order perturbation theory with 
T/Ec = 0.01,0.05,0.1,0.2,0.3,0.4,0.5,0.75,1,2,5, and 10. In this limit the scaled conduc- 
tance is independent of at- Right hand side: result of resonant tunneling with at = 0.2 and 
T/Ec = 0.01,0.05, and 0.1. 
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FIG. 3. The maximum linear differential conductance normalized to the high temperature limit 
for at = 0.001, 0.01, 0.1, 0.2 (from top to bottom). For comparison we also show the result obtained 
from lowest order perturbation theory (dashed line and inset). 




FIG. 4. The normalized nonlinear differential conductance for at = 0.1 and 
eV/Ec = 0.05,0.2,0.5 at zero temperature. The inset shows the result from a master equation in 
lowest order perturbation theory for eV/Ec = 0.5. 



CI(R,+RJ/a,e 




FIG. 5. The current-voltage characteristics of a symmetric (i?L = Rr) SET transistor in the 
strong tunnehng regime defined by eq. (|4^) at dt = 2.5 and difi'erent temperatures T/dt-Ec = 0.01 
(a), 0.2 (b), and 1 (c). For each temperature maximum and minimum currents are plotted. Curves 
with oscillations demonstrate the gate modulation efi'ect. This effect is completely suppressed in 
the high temperature limit (curve (c)). 



0.8 



0.6- 



G/G. 




0.5 



T/E,^ 



FIG. 6. Maximum and minimum linear conductance of a SET transistor as a function of tem- 
perature obtained from the Langevin equation analysis (eq. (^)) for dt = 2. The intermediate 
curve shows the linear conductance averaged over all values of the gate charge. Inset: Conductance 
as a function of the gate charge for the same dt at different temperatures T/Eq = 0.15 (a), 0.5 
(b), 1 (c) and 2 (d). 
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G/G, 




T/E, 

FIG. 7. The same as in fig. 6 for at = 3. The curves (a), (b), (c) and (d) in the inset correspond 
respectively to T/Eq = 0.02, 0.1, 0.5 and 1. 
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